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1 Introduction 



In this paper we are going to prove certain relations concerning Yangian for so(N) which has been 
formulated in [I] and to consider thoroughly corresponding numerical R-matrix defining Yangian for 
so(N) ~ spin(iV). 

Let A be a Lie algebra and T a (a — 0,1,2,...) be representations of A in spaces V a . Consider 
operators R ab (u) £ ~End(V a ®Vb), where u denotes a spectral parameter. We say that A is the symmetry 
algebra of the operator R a b (it) if Vg £ A we have 

(T a (g) ® h + la ® T b (gfj R ab {u) = R ab (u) (r a (g) ®I b + I a ® T b (gfj , 

where I a and I b are unit operators in V a and V b , respectively. Consider a set of Yang-Baxter RRR- 
equations for the operators R ab (u): 

R ab (u-v)R bc (u)R ab (v) = R bc {v)R ab (u)R bc (u-v) £ End(V <g> V b ® V c ) (1.1) 

where the representation spaces V a , V b , V c are different in general situation. There is an efficient pro- 
cedure which enables us to construct nontrivial solutions of cubic Yang-Baxter equations (|l.ip starting 
with the known one. The procedure can be illustrated by the following sequence of specializations in 
the Yang-Baxter relations (II. 1[) : 

Vb®Vb®Vb -4 V a ®Vo®V -4 V* ®V r a®V'o -> K«>K«)K -4 Vt,®V r (g>V r -> V 6 <8>Vb<g>V r -> H®^®^ 
and corresponding sequence of solutions 

Ro,0 - > Ra,0 — > Ra,a ^ Rfa.a ~> R b ,h ~> ' ' ' 

Indeed, one starts with the simplest known solution Ro.o of the Yang-Baxter equation (|1.1[) defined 
in the space Vo <8> Vo <8> Vb, where Vb is the space of the simplest faithful representation, e.g. defining 
representation for the matrix Lie algebra A. Further one introduces another representation T a which 
acts in the space V a (finite-dimensional or infinite-dimensional) and solves the Yang-Baxter equation 
(|l.ip restricted to V a €5 Vb <£> Vq. It happens to be a quadratic equation on the operator R aj o which in 
special cases represents the Yangian of the corresponding matrix Lie algebra A (Vb is the space of the 
defining representation and V a is the space of the representation of the Yangian) . On the next step one 
solves Yang-Baxter relation (jl.ip restricted to the space V a (8) V a ® Vq and obtains R a!<1 . There is a well 
known argumentation (based on the associativity ideas) why R a>a respects Yang-Baxter equation (|1.1[) 
defined in the space V a <8> V a <8> V a - Nevertheless it can be proved directly. Thus, the solution R a a of the 
cubic Yang-Baxter equation is constructed in several steps starting with the simplest one Ro,0: and in 
each step linear or quadratic relations have to be solved. 

To be more concrete let us remind [HE] now it works for the algebra A ~ so(d) ~ spin (cf), where we 
assume d to be even. All the following formulae can be rewritten straightforwardly for so(p, q) as well 
where p + q = d. Corresponding fundamental R-matrix R° (it) defined in the tensor product Vb <8> Vb of 
two fundamental (defining) (i-dimensional representations of so(d) can be represented as 

(R 0)H| (u)=u ^^ + ^^__|_^^ i . 2 ( (L2) 

and depicted as follows 

R°(u) 

R° respects Yang-Baxter equation 

R%(u-v)R° 3 (u)Rl 2 (v) = R° 3 (v)R° 12 (u)R° 3 (u-v) £ End(V ® V ® V ) (1.4) 



V, 



u + 



in 



(1.3) 
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and can be considered as the simplest solution in the hierarchy of solutions of the universal Yang-Baxter 
equation (jl.ip related to so(d) Lie algebra. It has been introduced in 1978 by A. Zamolodchikov and 
Al. Zamolodchikov [3J. 

On the next step we introduce spinor representation of so(d) acting in the space V with dimension 
2=. Let 7 (a = 1, . . . , d) be 2 5 -dimensional gamma-matrices in M d which act in V as linear operators. 
Operators j a represent generators of the Clifford algebra 

lalb +76 7a = 2£ a6 ' 1 ■ (1.5) 

As a vector space, the Clifford algebra has dimension 2 d . The standard basis in this space is formed by 
antisymmetrized products of the 7-matrices 

7a 1 ...a fc =^(-l) P(s) 7,(a 1 )---7,(a fc )^7A fc (Vfc < d) , lAk = (Vfc > d) , (1.6) 

where the summation is taken over all permutations s of fc indices {ai, . . . , afc} — » {s(ai), . . . , s(afc)} 
and p(s) denote the parity of the permutation s, is a multi-index a\ ...a/-. 

Then, according to procedure outlined above, we look for the operator L°(u) defined in the space 
Vo ® V which respects quadratic relation 

R° 23 (u-v)L 12 (u)L%(v)^L 12 (v)L 13 (u)R° 23 (u-v) G End(V ® Vb ® V ) , 

where R^u) is fundamental R- matrix (|1.2p . The solution of the above equation has been found in [2] 
(see also [1HS])- It has the form 

L»= U l®/„-i[7 Q ,7 & ]® e a& (1.7) 

where e a b are matrix units, 1 and /„ are identity operators in spinor and defining representation spaces 
respectively, summation over repeated indices is implied. 

Further from the universal Yang-Baxter equation (jl.ljl we obtain a linear equation for R-matrix 
Ri2(u) acting in the tensor product V ® V of two spinor representations 

Ri 2 (u- w)L? 3 (m)L° 3 (w) =L° 13 (v)L 23 (u)R 12 (u-v) g End(V ® V ® Vb) • (1.8) 

In [2] (see also 0HE]) spinorial R-matrix has been sought for in S'0((i)-invariant form 

R(u) = J-f^ 7a ^ ofc(87 a 1 ...a* G End(F(8V). (1.9) 

For convenience, in the r.h.s. of (|1.9p . the summation over k runs up to infinity. However we note that 
this summation is automatically truncated due to condition k < d (see (|1.6l0 . It has been claimed in [5] 
that the R-matrix (jl.9p satisfies RLL-relation (|1.8p if coefficient functions Rfc(w) obey the recurrent 
relation 

k — (u + d — 2) 

As far as we know it has not been checked directly so far that R(it) satisfies Yang-Baxter relation defined 
in the space V ®V ®V 

Ri2(w)R 23 (w + v)Ri2(v) =R 23 («)Ri2(u + w)R23(m) G End(V <g> V <g> V) (1.11) 

owing to complicated gamma-matrix structure one has to deal with. One of the aims of this paper 
is to carry out corresponding calculation. In order to avoid multiple summation over repeated indices 
we apply the generating functions technique and rewrite the sum in (|1.9j) as an integral over auxiliary 
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parameter. It enables us to perform the calculation in a concise manner. We undertake this calculation 
in Section [3] The derivation of the recurrent relations (jl.lOp is given in Appendix. 

Now we proceed to explore thoroughly spinorial R- matrix (ll.9[) . At first we are going to deduce its 
basic properties which can be obtained on a rather general argumentation and do not need intricate 
calculation technique. Further we formulate the other properties which we prove in subsequent Sections. 

It is well known that spinor representation T for generators M a b of the Lie algebra so(d) can be 
constructed out of gamma-matrices T(M a b) — §7ab f) 1 -6[) . Then one can easily check that R-matrix 
(|1.9j) is invariant under spin(d) action, i.e. 

(7 Qb <g> 1 + 1 <g> 7 Qb ) R(u) = R(u) {"f ab 1 + 1® 7a& ) , Va, b , 

moreover it satisfies commutation relation 

(7rf+i ® 7d+i) R(«) = R(«) (7d+i ® 7d+i) , (1-12) 
which demonstrates an additional u(l) symmetry of R(it). In (|1.12j) matrix 7d+i is defined as follows 

7d+i = «7i-d i « 2 = (-)* ; 7d+i = 1 ; {7rf+i:7a} = at a = l,---,d (1.13) 

and in the appropriate representation of gamma-matrices it takes the form 7^+1 = diag(/, —I). 

Let us note that recurrence equations (|1.10[) for the series of even coefficient functions R2k(u) and 
odd ones R,2k+i( u ) are independent. The general solutions of these equations are 

R 2fc (,)=AH(-l) fc r(fc + f^" fc) , M«)=B(«)(-1)*^|S^ (1.14) 



r(f)r(|) ' + v ; v JK ' 2r(^±i)r( 



where A and B are arbitrary functions of spectral parameter u. For example, if A and B are polynomials 
of spectral parameter then coefficient functions in (|1 . 14[) are normalized to be polynomials as well. Thus 
it is convenient to decompose spinorial R-matrix (|1.9I) in the sum R(u) = R + (u) + R~ (u) where (|1.6[) 

R+ («)=E^^®^ , R -(«)=E^§'TW®^ • (1-15) 

We refer to R + (u) and R _ (u) as even and odd parts of spinorial R-matrix, respectively. 
Consider the decomposition of the spinorial R-matrix in the sum 

R(u) = P+ R(u) + P" R(it) , (1.16) 

where are projectors 

P± = i (1 ® 1 ± 7d+1 ® 7d+1 ) ; P+P- = P-P+ = , (P±) 2 = P ± . (1.17) 

Proposition 1. Even and odd parts (| 1 . 1 5|) of the spinorial R-matrix can be singled out by projectors 
P^, i.e. we have 

R+(u) = P+R(u) , R-(u) =P-R(u) , (1.18) 
P ± R ± (m) = R ± (m) , P ± R =f (m)=0, R ± (u)R =F (w) = . (1.19) 
Proof. Due to (|1.14|) we see that coefficient functions in (|1.15|) satisfy the reciprocal conditions 

R 2fe (w) - (-l)*Rd- 2 fc(u) , Ra*+i(u) = -(-l)*Bd_ 2 fc-i(«) • (1-20) 
Taking into account (|1.6[) we deduce 

1 f 1) 2 

(7d+i ® 7d+i) 7j 7A fc ® l Ak = L lA d _ k ® 7^"* , (1-21) 
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(1.22) 



where Ad-k is the multi-index such that Ad-k H Ak = and A^-A; U At = {1,2, ... ,d}. Then using 
(fT20|) and ([TT2T]) we immediately obtain (|TTT8)) . Equations ([TTT9|) follow from ([LPT)) and (fTTHf . □ 

Proposition 2. The Yang-Baxter equation (jl.lip is equivalent to the following relations for R + and 
R" (fTT5l) : 

^•23 ^-12 = -R-12 -R-23 ^12 ! ^23 ^12 ^"23 = ^12 ^"23 ^12 ! 

R-23 R-12 -R-23 = ^12 -^23 ^12 ! -^23 ^12 ^23 = ^12 -^23 ^12 ! 
^23 ^12 ^23 = 7 R-23 ^12 -^23 = ^ ' ^23 ^12 ^23 = ; R-23 ^12 -^23 = ^ ' 
-R-12 ^-23 ^12 = , R12 R23 -R-12 = 1 R12 R23 -^12 = ; R12 P-23 ^12 = ) 

where the dependence on the spectral parameters is the same as in (jl.lip . 

Proof. The Yang-Baxter equation R^ 3 R+ 2 R23 = R12 R-23 ^12 i s deduced from (jl.lip if we act on it by 
projectors P^~ 2 and P^~ 3 from the left and right and use commutation relations 

p+ p + — p + p+ p+ T>+ — t>+ p + 
r 12 ^23 ~~ ^3 r 12 ! *23 rt 12 ~~ rt 12 r 23 • 

The relation R^ 3 R^ 2 R^ 3 = is obtained as following 

R-23 ^12 ^23 = ^23 ^12 ^23 ^23 = ^23 ^23 ^12 ^23 = ® ' 

where we use Rf 2 P^ 3 = P 23 R^ 2 , etc. □ 

We stress that in view of (jl.22l) the Yang-Baxter equation (jl.lip is satisfied for any linear combination 
R(u) = a(u)R + (it) + /3(w)R _ («) with arbitrary coefficient functions ot(u) and /3(u). It means that A(u) 
and B(u) in (|1.14|) are not fixed by equation (jl.lip . Moreover one can check by using (jl.22p that 
the Yang-Baxter equation (jl.lip is satisfied if we transform the solution R as following (we write this 
transformation in terms of even and odd parts R + (u), R~(u)): 

R+ -> R+ , R- -> ±R" ( 7d+1 ® 1) = T(7d+i ® 1) R- , 

(1.23) 

R+ -> R+ , R- -> ±R- (1 ® 7d+1 ) = t(1 ® 7rf+i) R • 
Proposition 3. Even and odd parts (j!.15p of the spinorial R-matrix satisfy unitarity relations 

R+(u)R+(-u) =h+(u)P+ , R _ (u)R~(-u) = h_(u)P~ . (1.24) 
where functions h + (u) , h_ (u) are constructed out of coefficents R^ (u) (|1.14p 

d/2 — — 1 

h+(u) = 2^ ( 2 d jR 2t («)R a (-«) = A(«)A(-«) I](fc 2 - u 2 ) , 

fc=0 ^ ' k=0 

d/2-1 , . |-1 

h_( U )=2 £ f jR 2fe+1 (u)R 2fc+1 (- U )=B( U )B(- M ) n( fc2 -" 2 )- 

fc=0 ^ ' fc=l 

Let us draw attention that in the right hand sides of the relations (| 1 . 24[) projectors P (| 1 . 1 7|) appear. 
Proof. At first in view of (| 1 . 14[) one obtains that at special value of spectral parameter spinorial R- 
matrix reduces to projector (j!.17|) : R + (e) = eT (|) P+ + 0(e 2 ) at e 0. Then the first Yang-Baxter 
relation in (jl.22|) at v = — u+e and e — > leads to Rj 3 (w) PJ2 R-23(~ u ) = ^i2(~ u ) R-23 R-i2(' u )- The latter 
relation is equivalent to P^ R 2 f 3 (w) R 2 f 3 (— w) = R^(— u) Rj 2 (w) R^ that implies R 12 (w) R^ 2 (— u) ~ P^. 
In a similar manner the second Yang-Baxter relation in (|1.22p leads to R 12 (u) R 12 (— u) ~ Pf 2 . Coefficient 
functions h + (u), h_(u) (jl.24p are calculated in Subsection [373] using generating function technique. □ 

Our considerations are aimed to the check of the Yang-Baxter equation (jl.lip for the R-matrices 
(jl.9p . (I1.14p and verification of their properties. For this we need to perform a rather complicated 
computations with Clifford algebra of gamma-matrices. To succeed in it we appeal to the technique of 
the generating functions developed in [7]. We briefly describe this technique in the next Section. 
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2 Clifford algebra 

2.1 Fermionic interpretation of Clifford algebra 

Let r a , a = 1, . . . , d, be a set of d generators of the Clifford algebra satisfying the standard relations 

(cf. <E3) 

r a r 6 + r 6 r a = 2<ui. (2.1) 

The Clifford algebra is a vector space with dimension 2 d . The standard basis in this space is formed by 
anti-symmetrized products of r a . For example (cf. (|f .6jl ) 

r.4o — 1 : r^i = r Q . Ta 2 = r aiQ2 = — [r ai r a2 — r a2 r ai ] , ■ • • (2-2) 



r Ak = r ai ... afc = As(r ai • • • r a j = A ^(-iy^r s{ai) ' ' • r ^ 



Here we use the notion of antisymmetric product As of r a -operators. Inside the As-product the operators 
T a behave like ant i- commuting variables. 

Note that the Clifford algebra generators can be represented as [IT] 

r a = 6 a + d 9a ; {9 a ,9 b } = 0, {d 8a ,d 9b } = 0, {dg a ,9 b } = 8 ab (2.3) 

where dg a — ^f- and 9 a (a = 1, • • • , d) form a set of d fermionic variables (generators of the Grafi- 
mann algebra). Below the fermionic interpretation of the operators T a will be important for us and to 
distinguish them from the matrices 7 a we use different notation ^ a — > T a . 
Now we introduce the generating function for the basis elements TA k (|2.2|) 



oo OO 

J2 £j u a " ■ ■ ■ « ai As (r ai . . . r a J = ^ - (u a r a ) k = exp ( u • r) = As [ exp ( u • r)] . (2.4) 

fc=0 ' k=0 

Here u ■ T = u a T a , u a are ant i- commuting auxiliary variables: u a u h = —u b u a and we also adopt that 
u a T b = —T b u a . Formula (12 .4[) implies that the basis elements r^ fe (|2.2p can be obtained from exp(u-r) 

as 

rai ...a k = d g exp( U -r) . (2.5) 

1 * u=0 

Further we indicate two basic relations which will be used extensively in our calculations with Clifford 
algebra. 

Proposition 3. The product of generating functions (12.41) is evaluated as 



e 



"i- r . . . e 1 ""' 1 " = e~ ^ ! <j Ui ' Uj e^ =1 U ^' F . (2.6) 



Let u a , v a , a a , f3 a be anti-commuting variables, x and y are commuting variables. Then we have the 
following identity 

exp (x d u ■ d v ) exp(u • a + v ■ (3 + y u ■ v) = (1 — xy) d exp ( a ■ (3 J , (2-7) 

u=u=o \1 — xy ) 



where we have used shorthand notation d u ■ d v = . 



Proof. The formula (|2.6[) is a consequence of the Backer-Hausdorff formula e A e B = e A + B +U A < B \ ) 
where we take A = u ■ T , B = v ■ T and [A, B] = -u a v b {T a T b + r b r a ) = -2u-v. 

Identity (|2.7|) can be easily deduced by taking into account the standard representation [8HTT] of the 
operator exp {xd u ■ d v ) as gaussian integral over 2d anti-commuting variables 9 a and 9 a . Indeed 

/d 
Y[ d9 a d9 a exp (a;- 1 9 ■ 9 + 9 ■ d u + d v ■ 9) , 
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so that all operations of differentiations lead to the simple shifts u — > u + 8 ,v — > v — 6 and then the left 
hand side of (|2.7I) takes the form of the gaussian integral again 

x d J Yl dO a d9 a exp ((x' 1 ~y)6-9 + 6-a + l3-9)=x d (aT 1 - y) d exp ^ - X ^ a ■ . □ 

In fact all subsequent calculations are based on (|2.6[) and (|2.7p . 

Note that the topic of this section has an evident interpretation in the language of quantum field 
theory. The formula (|2.6p is one of the variants of Wick's theorem and expresses the result of reduction 
to the normal form. The topic of this section can be considered as an application of the general 
field-theoretical functional technique [9 to a very special example, and exactly this point of view was 
elaborated in the paper [7 . It is possible to use the language of symbols of fermionic operators [TU] as 
well. For simplicity we have derived all needed formulae in a very naive and straightforward way. 



2.2 Fermionic realization of R-matrix 

Dealing with the Yang-Baxter equation (jl.lll) as well as with RLL-relation (|1.8I) we have to handle the 
tensor product of several spinor representation spaces. In fact we need gamma-matrices 7 a acting in 
the tensor product of two spaces. Since we consider instead of gamma-matrices 7 a the generators of 
Clifford algebra we need here two types of generators (Ti) a , (r 2 ) a which anticommute to each other 

(ri) B (r 2 )6 = -(r 2 ) 5 (ri) . (2.8) 

It is rather natural due to emphasized above fermionic nature of representation (|2.3p . Moreover the 
convention (|2.8I) makes the formulae much simpler. 

5'0(d)-invariant fermionic R-matrix (11.91) is constructed out of tensor products (Ti) Ak (r 2 )' 4fc . Let 
us rewrite this gamma-matrix structure in a more appropriate form 

(Ti) Ak (r 2 ) Afc = As [r lQl • • • r lak ] As [IT • • • rj*] = As [r lai • • • r lQ j r? ■ • • rj* = 

= A S(1) [r lai • • • r lafc • • • r?} = Sk a S(1) [(r x • r 2 ) fc ] (2.9) 

where Sk = (— l) - ^ ~ and we denote by As(i) the operation As applied only for the product of (Ti) a . 
Below we will omit index (1) in the notation As^) since for the expressions of the type (12.9[) we have 
As(i) = As( 2 ). At the first step in (|2.9[) taking into account definition (jl.6[> one can forget about one 
of the symbols As due to convolution of two antisymmetric tensors. Next it is possible to accomplish 
rearrangements taking into account that Ti ■ T 2 = — T 2 • IV The last equality in (|2.9p implies that 
As [e xri ' r2 ] is a generating function for the set of tensor products (ri)^ (T 2 )' 4fe 

Thus we have succeeded in rewriting the multiple summation over repeated indices in a compact form. 
Proposition. 

Asr e xrvr2 l = e *a..a. t^-H-r*! (2.11) 

L J In— v— v ' 

Proof. Using f|2.5p we obtain 
s fc Fi) Ah (T 2 ) Ak = s k a Uai • • • flL e" Tl d v a, ■ • • d v « k e v - r > = (d u ■ d v ) k e u ^ . (2.12) 



u— v— 



u—v—0 



Substitution of (|2~T21> into (|2~TU|) gives (|2~TT1) . □ 

Consider a fermionic analog of the operator (jl.9[> where coefficient functions are assumed to be 
arbitrary. Using generating function (|2.10p we represent this operator in several equivalent forms 

R(u) = g^ ( r l)Afe (r 2 )^=£M^^ Aa(e- r ^)| A!=0 = R(«|*)*Afl(e* r ^) , (2.13) 

k=0 ' fe=0 
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where we have used shorthand notation R(x) * F(x) = R(d x ) F(x)| . Note that all information about 
coefficient functions of the operator R in (|2.13l) is encoded in just one function R(u|x) 

R( MW ^^A (2.14) 

fc=0 

At the end of this Subsection we show how to represent fermionic operators (12.131) in the matrix form. 
There are two matrix representations p' and p" for the fermionic Clifford algebra with generators Y\ a 
and r 2a and defining relations (|2.ip . (|2.8[) : 



p'(Tla) = 7o ® 1 , P'(r 2 a) = 7d+l ® 7a 

(2.15) 

p"(Tla) = 7a ® Jd+1 , P"(r 2a ) = 1 ® 7a 



where j a , jd+i are standard 7-matrices defined in (jl.5[> and (| 1 . 1 3[) . For the even and odd parts of (|2.13|) 

R+ = E ^ (TO^ (r 2 ) A2fc , R" = E (^Tl)T CTx)^ (r 2 ) A2fc+1 , 

/c — k — 

we obtain by using (|2.15p the following representations 



(2.16) 



/ 3 '( R+ ) = E^fyj^ 2fc ®7 A2fc , p''(R-) = -^E(^^^ +1 ®^ 2fc+1 J( 1 ®^+i)- ( 2 - 17 ) 

Taking into account the fact that the solutions of the Yang-Baxter equation (jl.lll) admit transformations 
(jl.23[) we can use the following convention to construct matrix representation p of the Yang-Baxter 
solutions (|2.13p : 

^£^(^(^=£^^7^- (2.18) 



k\ v ,Ak v ' ) ^ k 

sk=0 / k=0 



Let us note that at even d 



pi (^As (Ti • T 2 )^ = p" (J^Ab (Ti • r 2 ) d ^ = 7 d+i ® 7d+i • (2.19) 
2.3 Exchange operators 

In this Subsection we examine simple examples of the operators presented in the form (|2.13j) . 
Let us consider the exchange operators P, P' defined by means of relations 

(r 2 ) a p - p (ro a ; (r^p' = P'(r 2 ) a . (2.20) 

We are going to show that it can be represented in the form (|2.13j) 

P = e x * As (e* riT2 ) - As (e ri ' r2 ) ; P' = e~ x * As (e xT ^) = As ( e - r ^) . (2.21) 
To be more concrete let us rewrite previous expression for operator P in the following form 



DO 



P = E i P2f k =£ E T ^ r ^ ■ ■ ■ ^ T 2^2 k - 1 " " F 

k=Q ' k=0 \a 1 <a 2 < -<a fc 
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The proof presented below will serve as a simple example to demonstrate typical calculations with the 
generating functions. Firstly we prove identities 

e s-r lAs ^ x r 1 .r 2 j =As L x r 1 .r 2 +s.{T 1+x r 2 )^ . As ( e *rvr 3 ) e t-r 2 = Ag J e xr I .r J +t.(r s+ «r 1 )j _ (2 22) 

The proof is rather simple and we perform it in detail for the first product. 

e^ 1 As (e xr ^) = e xd ^ e^ 1 e^i-H^ I 

V / \u— v— 

_ e xd u -d v e u-s+u-r 1 +vT 2 +s-r 1 i = As ( e x(r 1 +s)-r 2 +s-r 1 \ 

I u—v—Q \ J 

Here we apply successively (|2.1ip and (|2.6p . Thus (|2.22p is proven. In fact this calculation set the 
pattern for subsequent manipulations with generating functions. 

We rewrite equation (|2.20[) with the help of generating functions (|2.5I) . (|2.10[) 

P(x) * d Sa e s Tl As (e xivr2 ) | s=Q = P(sb) * d ta As (e xivr2 ) e tT2 | t=Q . (2.23) 

Substituting (12.221) in (I2.23P and calculating the derivatives with respect to s a and t a we obtain 

P(x) * As [(T u + xT 2a ) e Kri ' r2 ] = P(ar) * As [(T 2a + xT la ) e xT ^] , 

or equivalently 

[P(x) - d x P(x)} * As (r lQ e xV ^) = [P(x) - d x P(x)} * As (r 2a e*™) , (2.24) 

where in the last transformation we use the formula P(x) * x n F(x) = d T x l P(x) * F(x). As an evident 
consequence of (|2.24l) we obtain differential equation on the function P(x) 

d x P{x)=P(x) => P{x) = e x , 

that finishes the proof of (|2.21j) . 

In (|2.22[) we have found the product of the generating functions (|2.4p and (|2.10l) . It is exactly what 
we need to check RLL- relations (|1.8p . (|A.1|) giving rise to condition (|1.10p . Corresponding calculation 
is implemented in detail in Appendix. 

2.4 Generating function for Yang-Baxter and unitarity relation 

In this Subsection we examine thoroughly the gamma-matrix structure of the Yang-Baxter (jl.lip and 
unitarity (jl.24l) relations. In order to apply technique outlined above and in view of matrix represen- 
tation p (|2.18p we consider instead their fermionic analogues, i.e. the relations for fermionic operators 

We start with fermionic version of the Yang Baxter relation (jl.lip whose right hand side is a sum 
of operator tensor products 

(r 2 ) Afe (r 3 ) A » (ro^ (r 2 ) Sfe (r 2 ) Cfc (r 3 ) Cfe (2.25) 

multiplied by appropriate coefficient functions of spectral parameters. According to our approach instead 
of simplifying products of fermionic generators of Clifford algebra in (I2.25P we multiply corresponding 
generating functions (|2.10p depending on parameters x, y and z 

k^{e xV ^)k^{e v ^)^{ey v ^)^{\~xy) d k^e i ^ TM . (2.26) 

Expanding the latter formula into a series over x, y, z and picking out appropriate term one obtains 
()2.25p . Let us outline derivation of (|2.26p . Using (|2.1ip one can rewrite the product of the three 
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generating functions in (|2T2Uj) as follows e l9 »' a « e vd >' dt e zd v d i e" T2+ " T3 e s ri+ * r2 e p r2+ « r3 . Then due 
to (l2~6l) 

e u-r 2 +n-r 3 _ e s-ri+t-r 2 . e p-r 2 +g-r 3 _ e s-r 1 +(u+t+p)-r 2 +(v+q)-r 3 +t-u+p-t+q-v+p-u 

and applying several times (|2.7|) one obtains the desired result (|2.26j) . In much the same way generating 
function of the tensor product structure in the left hand side of the fcrmionic Yang-Baxter relation 
([i~TT|) has the form 

As (e^ 2 ) As (e^ 2 ^) As (e* r ^) = (1 - xy) d As (e&S r ^+^Sf ^r 3 +^ Fl . r ^ (2 2?) 

Let us mention that expressions (|2.26[) and (|2.27|) are almost identical. 

Dealing with unitarity relation (1 1 . 24|) for spinorial R-matrix we calculate RiafwlRiaf— u ) that forces 
us to consider tensor products of fermionic generators (J"i) Ak (T2) k {Fi) B (T2) k ■ Corresponding 
generating function is the following 

As(e xr ^)As(e^ r2 ) = (l-xy) d As(e^ r ^) . (2.28) 

Thus we have indicated the generating functions for tensor product structure of the relevant fermionic 
relations for spinorial R-matrix (j2.13[) . In the subsequent Sections using obtained results we will prove 
these relations. 

Remark. Equations (I2.26|) . (|2.27l) . (|2.28[) give the identities for the exchange operators (|2.21[) : 
P12 P23 P12 = P23 P12 P23 j P'12 P23 P'12 — P23 P'12 P23 ' 

pp = ^As(r 1 -r 2 ) rf , p'p' = t^-As(ri -r 2 ) d , pp' = p'p = 2 d i. 
2.5 Local Yang-Baxter relation 

Let us note that due to (|2.26j) and (|2.27[) the following local Yang-Baxter relation takes place 
(1 -xy^As^ 1 ^ 2 ) As(e zIvr3 ) As (e* ri - ra ) = 
= (1 - x'y')- d As (e x ' r * r A As (e 2 ' ri r2 ) As (e^ ^ 

where parameters x, y, z and x' , y' , z' are related by equations 

x + V = z' (1 + x'y') z(l + xy) = x' + y' z (x - y) = z' (x' - y') ^ 
1 — xy 1 — x'y' 1 — xy 1 — x'y' ' 1 — xy 1 — x'y' 

The last relation in (|2.30p and the product of the first two relations in (|2.30[) show that the functions 

_ z(x-y) _ z (x + y)(l + xy) 

(1 - xy) ' (1 - xy) 2 

are invariant under the transformation x, y, z — > x' , y' , z' . Thus, the points (x, y, z) and (x', y', z') lie on 
the curve C a .b defined by the equations 

z(x-y) = b(l- xy) 

(2.31) 

, (x + y)(l + xy) = a (x - y)(l - xy) 

where b = Ai and a = ^ are parameters which fix the curve. The geometrical picture is the following. 
The second equation in (|2.3ip defines the family of curves parameterized by a in the plane (x, y). Thus, 
it is possible to introduce new coordinates (x, y) — > (a, t) in the plane where t is a coordinate on the 



(2.29) 
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curve specified by a. The variable t is a coordinate on C a ,b as well. Then due to the first equation in 
(|2.3ip the coordinate z is determined by b and (x, y) or equivalently by b and (a, t). The transformation 
(x, y, z) — > (x' , y' , z') is equivalent to the change of coordinates t —> t' on the curve C a ^- Now we specify 
the coordinate t on the curve and chose, according to (|2.31j) . new variables (a,b,t) instead of (x, y, z): 

1 + xy x + y x — y x — y . „. 

a = , b = z , t = . (2.32) 

1 — xy x — y 1 — xy 1 + xy 

In terms of these new variables the transformation x,y,z x' , y\ z' looks very simple 

a — > a' = a , b — > b' — a , t t' = — . 

a t 

The t—¥t! transformation follows from the second relation in (|2.30[) which can be written as b/t — alt' . 

At the end of this Section we note that the local Yang-Baxter equations were introduced in [TJ] and 
applied to the investigations of 3d integrable systems in many papers (see e.g. [I31IT4] ). 

3 Yang-Baxter relation and unitarity 

In order to prove crucial properties of the spinorial R-matrix (|1.9|) we need to transform it to a more 
appropriate form. In Subsection 13.11 we rewrite spinorial R-matrix in fermionic realization (|2.13[) as an 
integral over auxiliary parameter 

R (") = r t^^* H«) As (e^) + b(u) As ( e -**vr a )] (3 .l) 
Jo (l + x^) " + 2 

where a(u) and b(u) are two arbitrary functions related with A(u) and B(u) appearing in ()1.14j) . Repre- 
sentation (13.11) happens to be very helpful since it enables to avoid multiple summations over repeated 
indices in (|1.9|) . Moreover the finite summation over k in (|1.9[) is substituted by an integral over auxiliary 
parameter. Thus the Yang-Baxter equation (|1.11[> which would assert equality of the two cumbersome 
multiple sums if we use representation (|1.9p . turns into an equality of two integrals. Using representa- 
tion (|3.1[) we check directly that the Yang-Baxter equation ([l.lljl is satisfied. More concretely we show 
that the equation is equivalent to the symmetry of a certain integral taken over the space of auxiliary 
parameters. 

3.1 Spinorial R-matrix 

Previously we have shown that gamma-matrix structure of spinorial R-matrix (|1.9[) can be simplified 
considerably using fermionic realization (|2.13p . Now we are going to make one more step rewriting 
the function K(u\x) in (|2.13[) that contains all information about coefficient functions Rfe(u). Let us 
remind that coefficient functions respect recurrence relations (jl.lOp . Above we have already found their 
solutions (I1.14p containing two arbitrary functions of spectral parameter. Using this freedom coefficient 
functions can be expressed in terms of Euler beta function 

R 2t M = A(-)'- r(t + ^ r( f - k) , .■.M^H'' 1 '^'^-' 1 (3.2) 

where A(u) and B(u) are arbitrary functions of spectral parameter. Then we separate even and odd 
terms in (12.14)) 



E R t(«) S k k ^ ^2k(u) S 2 k 2k , R 2fc+l(") S 2 k+1 2k+l 

k\ V ^ (2k)\ y ^ (2k + l)\ V 

k=0 fe=0 v ' fe=0 v ' 
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take into account s 2 k = S2fe+i = ( — ) k , resort to integral representation of the B-function 



r(fc + f)r(^ - fc) _ r°° dx x- 1 x 2k 
r(« + 1) ~ Jo (l + x 2 ) u +i 

r(k + «±i)r(s±^=i - fe) _ p dx x 11 - 1 x 2fe+1 



r(u + |) 7 (i + x 2 )«+^ 

and sum up the series obtaining 

/■°° (ir Irl" -1 

R ( u l^ = / m L + - ^ + B ) + - B ) e ~* V ] ■ ( 3 - 3 ) 

JO (l+Or)"' 1 " 2 

Thus we have managed to substitute finite set of coefficient functions appearing in (jl.9[) by the integral 
over auxiliary parameter. Finally, applying ()2.13|) we deduce the desired form (|3.ip of the spinorial 
R-matrix claimed above. In (|1.16[) we indicated natural decomposition of the spinorial R-matrix in 
the sum of even R + and odd R - parts (jl.151) . The formulae (|3.1[) and (J3T3J) imply the second natural 
decomposition 

R(u) = A(u)R+(u) + B(u)R-(u) = a(u) TZ + (u) + b(u) U~ (u) (3.4) 

where 

K+(u)^ f ^~ u l As(e^) , 1l-(u) = (°° **\*\^ As(e^) , (3.5) 
Jo (l+x^)" + 2 Jo (l + x^) n+ 2 

and a = A + B, b = A - B. 



3.2 Integral identity 



Now we are ready to establish the Yang-Baxter relation (11.111) . More exactly we will prove at first the 
Yang-Baxter relation for spinorial R-matrix in fermionic realization. Its tensor product structure has 
been already discussed in Subsection 12.41 To be more precise corresponding generating function for its 
right hand side (|2 .26[) and left hand side (|2 .27[) have been indicated. Then in the previous Subsection 
we have found out that coefficient functions of the spinorial R-matrix can be arranged in a sole function 
(I3.3|) . Further let us note that the Yang-Baxter relation (| 1 . 1 1 1) in fermionic realization is equivalent to 
the set of eight three-term relations for 1Z + , TlT (I3.5|) 

TZ\ 2 (u) K k 23 (u + v) n{ 2 (v) = TZ 23 (v) TZ k 2 (u + v) K 23 (u) (3.6) 

where k = +, — since a(u) and b(u) in the expression of spinorial R-matrix (|3.4p are arbitrary 
functions. At first the Yang-Baxter relation will be proven for 7l + (u). 

Taking into account (|2.26l) . (|2 . 27[) and (13.3[) one can easily see that the Yang-Baxter relation (|3.6[) 
at i = j = k = + is equivalent to 



As 



i"- , '(r 1 .r 2 ,r 1 -r 3 ,r 2 .r 3 ) 



As 



i"-"(r 2 -r 3 ,r 1 .r 3 ,r 1 .r 2 ) 



(3.7) 



where 



l u ' v (A,B,C) 



dxdydz \x\ u - 1 \y\°- 1 \z\"+ v - 1 (l-xy) d A ^ +B ^) +c£ M 

f 1 — my l — xy l — xy 



d (l+x 2 ) u +^(l + y 2 ) v +^(l + z 2 ) 



(3. 



the integration domain D = {(x, y,z):x>0,y>0,z> 0}. Instead of verifying (|3.7j) we are going to 
check more general relation 

F' V (A, B, C) = 1 U ' V {C, B, A) (3.9) 
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Figure 1: Projection of the domain D onto the plane (x, y). It is separated by the curves x = y, xy = 1 
into four parts marked by different colors, each mapped on the corresponding domain in Fig. [2j 

Figure 2: Projection of the domain G (|3.10|) onto the plane (t,a). It is separated into four parts 
corresponding to four subdomains on Fig. [TJ 



where the left and right hand sides to be understood as formal power series in A, B, C which are 
unspecified commuting external parameters. The discrete symmetry (|3.9I) of the integral (|3.8[) will be 
established by means of the integration variable change (x, y, z) — ► (x', y', z') defined by the local Yang- 
Baxter equation (|2.29|) which leads to the system of relations (|2.30[) . One can easily see that under this 
transformation of variables the external parameters A and C are interchanged in the exponential factor 
in (|3.8|) . However it is rather nontrivial that the other factors in the integrand transform in the right 
way such that Ij3.9|) is satisfied. 

To see it we appeal to geometric interpretation of the transformation (|2.30|) which have been discussed 
in Section l2"31 where we proposed to change variables (x, y, z) — > (a, b, t) according to (|2.32[) . In this case 
the integration domain D can be represented as (J fcCa.fc, where C a ,b is a curve parameterized by a and 
b. After all we make in (|3.8|) the natural change of integration variables (x, y, z) — > (a,b,t), presented 
in (|2.32l) . for which the Jacobian determinant has a rather simple form 



d(t, a, b) 



d(x,y,z) 



= 2 (l + x 2 )(l+ y 2 ) 



(1 + xy)(l - xy) 3 ' 

The formulae (|2.32[) map domain D onto disconnected domain G 

G = {(a, b, t) : a > 1, b > 0} U {(a, b,t) : a < -1, b < 0} (3.10) 

that is illustrated in Fig. [lj [2] After a simple calculation one obtains 

1 f dadbdt Ib^+v-^tl^+^+^aqyiAat- Bb + Cbt- 1 ) 
I U ' V (A,B,C) ^ ^ r / ^ J— i. 3.11 

2U+V ^ Jg |l + a\^\l - a\^- [6 2 + (1 + 6 2 )t 2 + a 2 f4]«+"+l 

Since the integral Q3.8P is rewritten in the form p. lip it is straightforward to prove the symmetry 
(|3.9j) applying the integration variable change (a,b,t) — > (a,b,t') where t' = It corresponds to 
the transposition at in p. lip . Indeed the integrand in p. lip transforms correctly and the 

integration domain G is mapped onto itself. 

Thus the Yang-Baxter relation p.6p at i = j = k = + is established. In a similar way the rest seven 
three-term relations p.6p can be checked. To realize it we note that expressions p.5p for 1Z + and 1Z~ 



12 



Figure 3: Projection of the domain D and of its three reflections onto the plane (x,y). Images of the 
subregions in the space (a, 6, t) are indicated (see (I2.32p ). It is assumed z > 0. If otherwise z < then 
b to be substituted on the figure by —b. 

are almost identical. They can be obtained from each other reflecting the integration variable x — > —x. 
In other words 1Z + and 7Z~ differ solely in integration contour. In the first case one integrates over 
positive semiaxis and in the second case over negative one. Consequently to check one of the three-term 
relations (13.61) we have to consider the integral (|3.8I) taken over appropriate reflected domain D. For 
example at i = j = — , k = + we integrate over x<0,y<0,z>0in (|3.8[) . When the variable change 
(|2.32p is performed, it leads to the integral p. lip with a certain integration domain that can be found in 
Figl21 The symmetry (I3.9[) is established as before by means of the variable change at <=^ bt^ 1 in p. lip 
which preserves the integration domain as one can easily see. Let us stress that algebraic manipulations 
needed to prove the three-term relations (|3.6p are the same in all eight cases. The only difference is in 
the integration domains in (I3.8P or (13. 

Finally, we have checked eight three-term relations (|3.6p and hence we have proved the Yang-Baxter 
relation (jl.lip for spinorial R- matrix in fermionic realization. Using decomposition (|3.4p of the R- matrix 
in the sum of even and odd parts we obtain eight three-term 

R[ 2 (u)R k 23 (u + v)R{ 2 (v) = R 23 (v)R k 2 (u + v)R 23 (u) (3.12) 

where k = +, — (compare with (| 1 .22[) V However let us emphasize that we have always used above 
the fermionic representation for the R-matrix. 

At the end of Subsection 12.21 we have shown how to represent fermionic operators (|2. 13|) in the 
matrix form (|2 . 18|) . It can be easily checked that three-term relations (|3.12p remain valid in both 
matrix representation p' and p" (|2 . 1 5|) . Thus the Yang-Baxter relation (jl.lll) for spinorial R-matrix 
(fL~9]) is checked. 

3.3 Unitarity relation 

In the Introduction we have formulated unitarity relations fll .24|) and proven them up to explicit cal- 
culation of coefficient functions h + (u), h_(u). Now we are going to fill this gap. We use fermionic 
realization of R-matrix. In view of (|2.13|) and (|2 . 28|) one has 

R + (tt)R+(-u) = R+ {u\x)R + \-u\y) * (1 - xy) d As (e^ ri ' r2 ) . (3.13) 

Since we know that R + (u)R + (— u) is proportional to projector P + , formula (|3.13p contains only fermionic 
structures 1 and As (r! ■ Y 2 ) d . Coefficients at the other structures are equal to zero. Thus it will be 
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sufficient for us to caicuiate numericai coefficient for 1 in p. 131) that is equai to 

R+(u\x)R+(-u\y) * (1 - xy) d = £^ 2 ( 2 d jR 2fe (w)R 2fc (-7i) . 

Similariy calculating numerical coefficient for As (ri • T2) d in (|3.13p one obtains 

R+(u|x)R+(-u|y) * (x + y) d = E^o(-l)* ( 2 d k )R2 k (u)R d - 2k (-u). 

Further, using matrix representations p 1 (|2 . 16[) or p" ()2.17j) one obtains (|2.18|) 

p'(R + (u))p'(R + (-u)) = p"(R + (u))p"(R + (-u)) = p(R+(«))p(R+(-u)) 

that leads finally to the first unitarity relation (ll.24j) in view of (|2.19[) . 

The previous arguments are valid also for R _ (u)R~(— u). The coefficient function for 1 is equal to 

R-(u\x)R-(-u\y) * (1 - xy) d = - E^Jq 1 ( 2k d +1 )^2k+i(u)R2 k+ i(-u) . 

The matrix realization of the second unitarity relation (|1.24|) is provided by 

p'(R-(u))p'(R-(-u)) = p"(R-(u))p"(R-(-u)) = -p(R-(u))p(R-(-u)) . 

Remark. Unitarity relations (|1.24[) can be established as well be means of integral representation for 
R-matrix (|3.1|) using 



oo oc 

// 



dxdyx u 1 y u 1 (x + y) k (I — xy) d k 



u sin 7TU 
2tt 



at k = 







cot iru , at k = d 
0, at k = 1, - ■ ■ ,d — 1 . 
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A Appendix 

In pQ we introduced representation space V' which we assumed to be infinite-dimensional in general 
and restricted the universal Yang-Baxter equation (jl.ip to the space V <S) V (g> V 

R 12 (u-v)L 13 (u)L 23 (v) =L 13 (v)L 23 (u)Ri 2 (u-v) e End(U ® V ® V) . (A.l) 

The operator L(u) which is defined in the tensor product V ®V of spinor and arbitrary representation 
T" spaces has been sought for in the form 

L(tt) = ti+ l - lab ®T'{M ab ). (A.2) 

Here notation (12.21) is used and M a b (a,b — 1, . . . , d) are generators of so(d) subjected to relations 

[M ab , M dc ] = i(S bd M ac + 5 ac M bd - 5 ad M bc - 8 bc M ad ) . (A.3) 
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In pQ we claimed that RLL-relation (|A. 1|) with spinorial R-matrix (|1.9p is satisfied if representation T" 
is such that 

T'({M [ab ,M c]d }) =0, (A.4) 

where {j4,J3} = AB + £?A is anticommutator and square brackets denote antisymmetrization. We 
will undertake corresponding calculation in the first part of this Appendix using generating function 
technique. We are going to show that RLL-relation (|A.1|) with L-operator (|A.2[) and R-matrix of the 
form (|1.9j) leads to recurrence relation (|1.10j) for coefficient functions R& (it) and set up restriction (|A.4|) 
on the representation T' in the quantum space. 

In order to avoid misunderstandings let us note that in a special case d = 6 we have the isomorphism 
so(6, C) = si(4,C), corresponding 8-dimensional L-operator (|A.2j) is a direct sum of two 4-dimensional 
L-operators of s£(4) algebra, spinorial R-matrix (ll.9[) reduces to Yang R-matrix under Weyl projections 
and condition (|A.4[) on representation T" happens to be superfluous. We demonstrate it in the second 
part of this Appendix. 

A.l RLL-relation 

Further by abuse of notation we denote T'(M ab ) — > M ab . The following calculation is very similar to the 
one presented in Subsection 12. 31 and it uses the generating function technique. We are going to prove 
fcrmionic version of RLL-relation (|A.1[) . Then taking matrix representation p' or p" (|2.15p one obtains 
immediately (|A.1[) for spinorial R-matrix (jl.9[) and L-operator (|A.2[) . 

The substitution of spinorial R-matrix (|2.13[) in fermionic realization and fermionic analogue of 
L-operator (|A.2j) with unspecified representation T' in the quantum space in RLL-relation (|A.1|) gives 



k\ 

k=0 



cd 



= E ?i ir( + i( r ^) (- + ^(r 2 ) cd M cd )(r 1 ) Afc (F 2 )^ . (A.5) 

This relation contains terms linear and quadratic in generators M ab . The product of two generators can 
be transformed by means of Lie algebra commutation relations (|A.3|) 

lli 1 

M ab Md = - [M ab , Med] + - {Mab , M cd } = - [g bc M ad - g ad M cb - g ac M bd + g bd M ca \ + - {M ab , M cd } 

so that 

(TiU (r 2 ) cd M ab M cd = -2i (r x ) a c (r 2 ) bc m ab + \ (r x ) a6 (r 2 ) cd {M ab , M cd } . 

All terms in (|A.5|) linear on spectral parameters are combined in a single one ~ (u — v) due to relation 

(Ti) Ah (T 2 ) A * (T 2 ) ab - (T,) ab (r x )^ (T 2 f- = (r x )^ (F 2 ) afc (T 2 ) A « - (T,)^ (T^ (T 2 ) A « , 
that is a consequence of the so(d) invariance 

(ri) o6 + (r 2 ) o6 ,<Ti) At (r 2 ) Ak 

After all intertwining relation (|A.5[) is reduced to the form 

E u M ° 6 (Vi W (r 2 ) Afc (r 2 ) ab - (T^ (T^ (F 2 )^ ) + 
+^ E ^# Ma " (Vi W ( r i) Q c ( r ^ fc ( T *) bc (ri) a c (rr) Afc (r 2 ) fcc (r 2 )^ ) + 
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i g Rjp_ | (r2)Afc ^ _ (ri) ^ ^ (r2)Afc | {Ma6 ^ Mcd} = Q _ (A _ g) 



fc=0 



Using the reference formulae for products of generating functions (see (|2.22[) ) 



As 



As 



As 



(e* ri 



r 2 +t-(r 2 +xr 1 ) 



s'^ 2 As 



it is easy to derive compact expression for the first term in (IA.6I) 



£ |m»* (ro^ (r 2 ) A * (r 2 ) ab - (r^ (r^ (r 2 y 



fc=0 



R{x) * M ab d Sa d Sb Ase xri r2 



e s-(r 2 +xri) _ e s-(ri+ar 2 ) 



R(x) * M^Ase^ 1 ^ 2 



(r 2a + ^ri b )(r 2fc + xT la ) - (T la + xr 2a )(ri h + xT 2 b) 



= R{x)*(x 2 - 1) M ab As e xTl ' T2 [T la T u - Y 2a Y 2b \ = (d 2 R{x) - R{x))*M ab Ase xT ^ 2 [T la T lb - T 2a T 2b ] 
In a similar way using 

As(e xivr2 ) e s ' ri+t ' r2 = As (e a:(ri_s) ' (r2_t)+s ' ri+t ' r2 ^ 
e s-ri+t-r 2 Ag (gzri-r 2 j _ Ag ^ e x(r 1 +s)-(r 2 +t)+s-r 1 +t-r 2 'j 
the second term in (|A.6[) can be rearranged as follows 

£ M ab ((ro^ (T!) a c (r 2 ) A * (r 2 ) bc - ( ri ) a c (r^ (r 2 ) hc (r 2 ) A *) = 



2RW * M Qb Ase* ri r2 [r la r lb - r 2Q r 2b ] 



(x 3 + x) r lc r c 2 - (d - 2) x 2 



= -2 



x dlR{x) + x d x R{x) -{(1-2) d 2 x R{x) 



M a b Ase *rvr 2 [r la r lh - r 2a r 2b ] 



and the last term in (|A.6[) takes the form 



E if ( (^W (ri) ab (r 2 ) Afc (r 2 ) cd - (r x ) cd (r x ) A (r 2 ) ah (r 2 )^ ) {m° 6 , M cd } 



fe=0 



4R(i) * (x 3 - x) {M ab , M cd ) Ase xri ' r2 



riaTibrici^d _ r 2a r 2b r 2c ri ( i 



9 3 R(x) - d x R(x) 



{M ab ,M cd }As e xTl ' T2 



riaribricr 2 d — r 2a r 2b r 2c rid 



Thus finally we obtain that (|A.6|) is equivalent to the relation 



xdlR{x) +xd x R(x) -(d-2) d^R(x) - u {d 2 x R{x) - R{x)) 



*M ab Ase xr ^ [T la T lb - T 2a T 2b ] 



d 3 R(x) - d x R(x) 



* {M ab ,M cd } Ase xIvr2 



riarihTici^d — r2ar2;>r2crid 



(A.7) 



l(i 



There are two independent gamma-matrix structures in the latter formula so that the differential equa- 
tion for the coefficient function R(x) 



fl£R(ar) + d x R(x) 



-(d-2)d%R(x) - u 



d*R"(x)-R{x) 



and requirement {-Mj afc ,M c ] d } = (see (|A.4I) arise. The differential equation produces the recurrence 
relation (see (|1.10p ) for the coefficients Rfc(tt): 



fe=0 



fc! 



R, 



fc+2 



(U) = - 



d-2-k 



R*(u) 



A. 2 R-matrix in the special case d = 6 

Now we proceed to the special case d = 6. The recurrent relations (|1.10[) for odd and even coefficients 
are independent that enables us to fix Ro(u) = (u + 4)/8 and Ri(u) — 0. Hence R-matrix (11.9[) takes 
the form 



R(u) = Ho(u) 1 ® 1 + 5^) 7aiQ2 ® 7 — + 7ai ... 04 g, 7 ax...«« + MO 7qi ae ? 

where 

Ro(it) = (u + 4)/8 , R 2 (it) = -it/8 , R 4 (w) = u/8 , R 6 (u) = -(u + 4)/8 
and the last term in (IA.7j) which is responsible for the condition (|A.4|) reduces to 

2 " 



7" 



(A.8) 



3! 



R 6 (u) +R 4 (u) 



{M o6 ,M° d } 



7a6cc 1 c 2 c 3 ^7 i / ^ ;afccc 1 c 2 C3 



(A.9) 



All the other terms vanish because of the special form of coefficients Rfe(u) and owing to finiteness of 
the Clifford algebra of gamma-matrices. Next we note that owing to ai abcciC2C3 = e abcciC2C3 7 7 and 
7777 = 1 (|1-13[) the gamma-matrix structure in (|A.9[) can be transformed as follows 

labcc lC2 c 3 ® 7 dClC2C3 = 77 <8> 77 7abcc lC2 c 3 7 dClC2C3 = 120 77 ® 77 [ $t Jbc ~ # lac + Sf lab } . 

Consequently (|A.9[) which is proportional to 



{M a \M c d } [5 d albc -5tiac + 5 d clab ] =2{M a ( b ,M c \} lbc 







turns to zero. In the last expression the parentheses (...) denote symmetrization. Therefore RLL- 
equation (|A.5[) is valid for arbitrary representation of generators {M ab } of the algebra so(6). 

Let us rewrite the expression for R-matrix (|A.8|) in a more transparent form. All gamma-matrix 
structures in (IA.8[) have block-diagonal form in Weyl representation for gamma-matrices. Therefore 
it is reasonable to consider projections of (|A.8[> on corresponding irreducible subspaces. We introduce 
subspaces V+ and V- obtained by Weyl projections: V+ = 1+ F 7 V and V- 



l-r, 



V . At first we note 



that relations 



1_ &. lAe 



n 



i 



l 



lead to R{u) 



V+®V- ~ 

1®1 



= R ( u )\ V -®v+ 

1 

+ gj 1A 6 ® 7 



- 7A 2 ® 7 - -77 7A 4 

0. Further a pair of relations 



v^4 



= 



V+0V- 



V-&V- 



1A 2 



leads to Yang R-matrix 

R(")ly_®y_ = [2R (u)l® 1 + R 2 (w)7afc®7 ab ] 



V_(8V_ 



where P is a permutation operator and we take into account — g 7a f, ® 7 ob | ^ y. 
Analogously one concludes that R(u)| v „ v =l®l + uP. 



= P-il 
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